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摘 要 : [目的 /意义 ] 人 工 智 能 (Artificial Intelligence, AI) 技术 已 在 学 术 和 工程 应 用 领域 掀起 了 研究 高 淹 ， 在 
地 球 物 理 参数 和 农业 气象 遥感 参数 反 演 方面 也 表现 出 了 强大 的 应 用 潜力 。 目 前 大 部 分 AI 技术 在 地 学 和 农学 的 应 用 
还 是 “黑箱 ”， 没 有 物理 意义 或 缺乏 可 解释 性 及 通用 性 。 为 了 促进 AI 在 地 学 和 农学 的 应 用 和 培养 交叉 学 科 的 人 才 ， 
本 研究 提出 基于 AI 耦合 物理 和 统计 方法 的 地 球 物理 参数 反 演 范式 理论 。[ 方 法 ] 首先 基于 物理 能 量 平衡 方程 进行 
物理 逻辑 推理 ， 从 理论 上 构造 反 演 方程 组 ， 然 后 基于 物理 推导 构建 泛 化 的 统计 方法 。 通 过 物理 模型 模拟 获得 物理 
方法 的 代表 性 解 以 及 利用 多 源 数据 获得 统计 方法 代表 性 的 解 作为 深度 学 习 的 训练 和 测试 数据 库 ， 最 后 利用 深度 学 
习 进 行 优化 求解 。[ 结 果 和 讨论 ] 判定 形成 具有 通用 性 和 物理 可 解释 的 范式 条 件 包括 : (1) 输入 与 输出 变量 CE 
数 ) 之 间 必 须 存在 因果 关系 ; (2) 输入 和 输出 变量 (参数) 之 间 理 论 上 可 以 构建 闭合 的 方程 组 (未知 数 个 数 少 于 
或 等 于 方程 组 个 数 )， 也 就 是 说 输出 参数 可 以 被 输入 参数 唯一 确定 。 如 果 输 入 参数 (变量) 和 输出 参数 (变量 ) 之 
间 存 在 很 强 的 因果 关系 ， 则 可 以 直接 使 用 深度 学 习 进 行 反 演 。 如 果 输 入 参数 和 输出 参数 之 间 存 在 弱 相 关 性 ， 则 需 
要 添加 先 验 知识 来 提高 输出 参数 的 反 演 精度 。 此 外 ， 本 研究 以 农业 气象 遥感 中 的 关键 参数 地 表 温 度 、 发 射 率 、 近 
地 表 空 气温 度 和 大 气 水 汽 含量 联合 反 演 作为 案例 对 理论 进行 了 证 明 ， 分 析 结 果 表 明 本 理论 是 可 行 的 ， 并且 可 以 辅 
助 优化 设计 卫星 传感器 波段 组 合 。[ 结 论 ] 本 理论 和 判定 条 件 的 提出 在 地 球 物理 参数 反 演 史上 具有 里 程 碑 意 义 。 
关键 词 : 人 工 智 能 ; 深度 学 习 ; 反 演 范式 ; 物理 逻辑 推导 ; 可 解释 ;农业 气象 遥感 
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1 8| &E 受到 人 们 的 广泛 关注 ，AI 成 为 新 一 轮 产 业 变革 的 核 
心 驱 动力 ， 将 进一步 释放 科技 创新 的 巨大 能 量 ， 并 

以 ChatGPT 为 代表 的 生成 式 人 工 智能 (Artifi- ”创造 新 的 强大 引擎 。AI 与 具体 的 行业 结合 将 催生 新 
cial Intelligence, AI) 的 出 现 和 升级 使 得 AI 进 一 步 。 技术 和 新 产品 ， 深 刻 改变 人 类 思维 和 生产 模式 ， 实 
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现 社会 和 行业 生产 力 的 整体 跃迁 。 深 度 学 习 是 目前 
AI 领域 最 重要 的 技术 之 一 ， 在 学 术 领 域 和 工程 应 用 
掀起 了 研究 高 潮 ， 并 在 自然 语言 生成 模型 、 计 算 机 
视觉 、 语 音 识别 等 多 个 领域 取得 了 显著 的 研究 成 
果 。 鉴 于 深度 学 习 在 地 学 和 农学 领域 的 应 用 潜力 和 
重要 性 ， 为 促进 AI 在 地 学 和 农业 领域 的 应 用 ， 必 
须 加 快 AI 与 遥感 技术 的 深度 融合 ， 为 气象 预报 、 
农 情 监 测 和 农业 灾害 预报 等 提供 关键 技术 支撑 S, 
为 全 球 灾害 监测 和 国家 粮食 安全 服务 。 

根据 2022 年 《国家 统计 年 鉴 》， 中 国 年 均 农 业 
气象 灾害 受灾 面积 达 1207.16 HAW, 122022 年 气 
象 灾害 造成 的 直接 经 济 损失 就 高 达 2386 亿 人 民 币 。 
因此 ， 快 速 准确 获取 农业 气象 关键 参数 信息 对 于 科 
学 指导 农业 气象 防 灾 减灾 、 作 物 估 产 确保 国家 粮食 
安全 等 具有 重要 意义 “*"。 美国 、 日 本 和 欧盟 等 相 
继 研 发 气象 卫星 关键 参数 反 演 系统 ， 中 国 自 20 志 
纪 90 年 代 开 始 研究 气象 卫星 参数 反 演 技术 ， 虽 然 
都 已 经 取得 了 很 大 的 进展 ， 但 仍然 遇 到 一 些 瓶 颈 ， 
如 未 知 数 多 于 方程 的 个 数 、 相 对 和 人 射 角 变化 和 混合 
像 元 等 问题 ”。 近 年 来 ， 在 计算 机 硬件 技术 和 AI 
技术 迅速 发 展 相互 辅 助 下 ，AI 显 示 出 了 强大 的 优 
势 ""。 因 此 研究 如 何 将 AI 技术 与 地 学 ， 特 别 是 与 


地 球 物理 参数 的 反 演 精度 。 这 样 保证 了 反 演 范式 不 
仅 具 有 通用 性 ， 而 且 具 有 物理 意义 和 可 解释 性 。 具 
体 地 ， 本 研究 提出 的 地 球 物 理 参数 反 演 范式 理论 是 
使 深度 学 习 的 输入 参数 (变量 ) 和 输出 参数 (AB 
量 ) 之 间 可 以 构建 一 套 完 整 的 闭合 方程 组 。 如 果 输 
入 参数 (Ait) 和 输出 参数 (变量 ) 之 间 存 在 强 因 
果 关 系 ， 则 可 以 直接 使 用 深度 学 习 进 行 反 演 并 获得 
高 精度 的 反 演 结果 。 如 果 输 入 参数 和 输出 参数 之 间 
存在 弱 相 关 性 ， 则 需要 添加 强 相关 性 的 先 验 知识 和 
交叉 迭代 来 提高 输出 参数 的 反 演 精度 。 可 以 通过 物 
理 模 型 模拟 得 到 物理 方法 的 大 量 代 表 性 解 ， 然 后 使 
用 深度 学 习 通 过 训练 获得 解 的 曲线 函数 。 因 此 通过 
大 数据 模式 的 深度 学 习 可 以 重复 该 物理 方法 。 同 时 
很 多 物理 方法 不 能 描述 所 有 情况 ， 为 了 克服 物理 方 
法 的 不 足 ， 基 于 物理 逻辑 推理 ， 进 一 步 构建 泛 化 的 
统计 方法 ， 并 通过 多 源 数据 获得 统计 方法 的 解 来 弥 
补 物理 方法 的 不 足 。 其 实 从 信息 流传 递 的 角度 来 
说 ， 深 度 学 习 网 络 更 能 符合 现实 情况 ， 物 理 方法 和 
统计 方法 是 其 中 的 部 分 表现 形式 。 这 些 方 法 本 质 上 
一 致 ， 均 是 人 们 用 来 认识 世界 的 一 种 手段 ， 在 不 同 
认 知 程度 上 解读 和 解决 现实 问题 。 

判定 基于 深度 学 习 的 地 球 物理 参数 反 演 形成 是 


农业 气象 还 感 参数 传统 反 演 物理 和 统计 方法 结合 ， 
发 展 具 有 物理 意义 和 可 解释 的 AI 算法 对 提升 遥感 
参数 反 演 精度 和 对 完善 地 球 物理 参数 反 演 范式 理论 
具有 重要 意义 ， 同 时 为 智慧 农业 和 智慧 地 球 提供 关 
键 的 监测 参数 。 


2 ”地球 物理 参数 反 演 范式 理论 与 判定 
条 件 


建立 学 科 体系 知识 一 般 遵 循 三 个 步骤 : 首先 从 
物理 角度 观察 现象 ， 其 次 用 数学 方法 描述 物理 现 
象 ， 最 后 从 工程 的 角度 去 解决 数学 问题 。“AI 地 球 
物理 参数 反 演 范式 理论 ”可 以 表述 为 : 如 果 目 标 信 
A (问题 ) 在 理论 或 物理 逻辑 上 可 以 用 数学 方程 描 
述 ( 即 可 以 在 空间 中 可 以 形成 唯一 的 解 曲线 )， 那 
么 AI 方法 可 以 通过 大 数据 学 习 和 优化 计算 耦合 物 
理 方法 和 统计 方法 求解 方程 。 这 使 得 耦合 方法 不 但 
具有 物理 方法 和 统计 方法 各 上 自 的 优势 ， 同 时 充分 利 


否 具 有 通用 性 和 物理 可 解释 范式 的 条 件 有 两 个 : 

(1) 深度 学 习 的 输入 与 输出 变量 (参数 ) 之 间 
必须 存在 因果 关系 ， 即 在 其 他 条 件 不 变 的 情况 下 ， 
改变 任何 一 个 输入 变量 ， 输 出 变量 会 发 生变 化 。 设 
深度 学 习 输 入 变量 为 x， 输 出 变量 为 yy， 则 因果 关 
系 表示 为 公式 (1)。 


y,7 f(x) (1) 
HB, i=l, 2, n; j=l, 2, +, m; [fR 
示 某 种 因果 关系 函数 。 


(2) 深度 学 习 的 输入 和 输出 变量 (参数 ) 之 间 
理论 上 可 以 构建 闭合 的 方程 组 (未知数 个 数 少 于 或 
等 于 方程 组 个 数 ) ， 也 就 是 说 输出 参数 可 以 被 输入 
参数 唯一 确定 ， 用 数学 方程 可 以 表示 为 公式 (2). 

F(x, y) =0 (2) 

其 中 ，x =x, Xp. °°, X3 VEY Vor Uo X6 
FRR TH eT Zr EB E t PR AL 

为 了 使 输出 参数 可 以 被 输入 参数 唯一 确定 ， 需 


用 了 深度 学 习 的 优化 计算 能 力 ， 最 大 限度 地 提高 了 


要 满足 E> m 这 个 条 件 。 为 论证 深度 学 习 的 使 用 具 
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有 物理 意义 ， 在 研究 时 必须 进行 物理 逻辑 推理 。 输 
入 变量 和 输出 变量 之 间 存 在 因果 关系 且 理 论 上 也 能 
构造 足够 的 方程 组 时 ， 如 果 输 入 变量 和 输出 变量 之 
间 存 在 强 相关 性 ， 则 可 以 进行 直接 反 演 计算 且 可 以 
保持 高 精度 ; 但 如 果 输 入 变量 和 输出 变量 之 间 是 弱 
关系 ， 要 保证 获得 较 高 精度 则 需要 利用 强 相 关 性 变 
量 作为 先 验 知识 。 

深度 学 习 的 应 用 不 能 脱离 具体 的 问题 或 对 象 ， 
如 图 1 所 示 ， 本 研究 在 解决 某 个 应 用 问题 时 ， 首 先 
要 对 问题 对 象 进行 物理 逻辑 推理 ， 从 理论 上 构造 方 
程 ， 确 定 输入 变量 和 给 出 变量 之 间 的 因果 关系 ; 其 
次 从 方程 里 面 的 未 知 数 个 数 确定 方程 组 的 个 数 ， 也 
就 是 确定 几 个 输入 参数 能 唯一 确定 输出 变量 ; 再 次 


通过 分 析 确 定 利用 物理 模型 是 否 可 以 模拟 得 到 所 有 
代表 性 的 解 ， 如 果 得 不 到 ， 则 构建 基于 物理 逻辑 推 
理 的 泛 化 的 统计 方法 “""”， 利 用 多 源 数据 补充 代表 
性 的 解 ; 然后 利用 大 数据 技术 和 深度 学 习 优化 求解 
物理 和 统计 方法 解 的 空间 曲线 函数 ; 最 后 进行 验证 
和 反复 和 迭代， 达到 精度 后 停止 优化 计算 ^07 dn 
果 应 用 深度 学 习 时 ,没有 证 明 或 不 能 证 明 存 在 上 述 
两 个 条 件 ， 则 一 般 认 定 深度 学 习 的 应 用 为 “ 黑 
箱 ”"“"”"。 如 果 能 证 明 输 入 和 输出 参数 之 间 肯 定 不 
存在 闭合 的 关系 未知数 个 数 大 于 方程 组 个 数 )， 
那么 这 样 的 应 用 只 适合 于 局 部 区 域 ， 不 能 通用 和 很 
难 具备 移植 性 ， 因 此 不 能 称 为 范式 。 


基于 人 工 智能 的 地 球 物理 参数 反 演 范式 理论 及 判定 条 件 


物理 方法 


B(T,) = &zr,B(T,)- (1 c) [I- 7, (1-85)]8,(7,.7,) (=4) 


i 
物理 模型 模拟 解 
CD) yix) MODTRAN smulation data 
FHI, 2, m. 1,2, ms /表示 函数 
(2) Fix, )-0 
ik Pex 


范式 理论 判定 与 约束 条 件 


Ft 


XN is X °°. Ret PEV Bae te Mad e 
示 一 个 由 上 个 方程 组 成 的 向 量 函数 ， 且 
>m 


深度 学 习 耦 合 物理 和 统计 方法 
(农业 气象 通 感 关键 参数 反 演 范式 案例 分 析 》 


"s T 


地 表 温度 、 发 射 率 、 近 地 表 空 气温 度 
大 气 水 汽 含量 等 
ET . zi 


实测 数据 


图 1 地 球 物理 参数 反 演 范式 理论 与 判定 条 件 


Fig. 1 The paradigm theory and judgment conditions of geophysical parameter retrieval 
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3 地球 物 理 参 数 反 演 范式 理论 案例 分 
析 论 证 


3.1 物理 逻辑 推理 


农业 气象 遥感 关键 参数 是 地 球 物理 参数 的 重要 
组 成 部 分 ， 本 研究 以 农业 气象 遥感 关键 参数 为 例 对 
地 球 物理 参数 反 演 范式 理论 进行 论证 。 

首先 进行 物理 逻辑 推理 。 地 表 热 辐射 在 通过 大 
气 达 到 卫星 传感器 的 过 程 中 ， 主 要 受 地 表 类 型 和 十 
壤 水 分 (Soil Mositure, SM) 、 地 表 温 度 (Land 
Surface Temperature, LST) 、 近 地 表 空 气温 度 
(Near Surface Air Temperature，NSAT) 和 大 气 水 汽 
含量 (Water Vapor Content，WVC) 等 的 影响 。 地 
表 温 度 和 土壤 水 分 反 演 是 基于 地 表 热 辐射 传导 及 其 
通过 大 气 到 达 传 感 带 的 传送 物理 过 程 ， 其 反 演 方程 
通常 可 以 用 公式 (3) 描述 ， 如 图 2 所 示 。 


注 : 参数 为 土壤 水 分 、 地 表 温 度 、 近 地 表 空 气温 度 和 大 气 水 汽 含量 
图 2 农业 气象 逐 感 参数 之 间 的 关系 
Fig.2 The relationship between agricultural meteorological re- 


mote sensing parameters 


公式 (3) P, B, (T). (已 知 数 ) 为 卫星 接收 
的 星 上 辐射 ，(1 (0))(1+(1-&)w(0))B,(7T,) 为 
大 气 辐射 贡献 ，B,(T,) (0)s; 为 地 表 辐 射 。 其 中 ， 
(0) 为 大 气 透 过 率 (未 知 数 )，T 为 地 表 温 度 (未 
知 数 )， 工 为 近 地 表 空气 温度 (未知数 )，s, 为 地 表 
发 射 率 (未 知 数 )， 即 一 个 方程 至 少 有 4 个 未 知 数 。 
需要 说 明 的 是 观测 角度 其 实 也 是 一 个 未 知 数 ， 每 个 
像 元 中 不 同 地 物 的 组 合 导致 相对 入 射 角 变 化 很 大 ， 
很 难 被 准确 地 确定 。 由 于 地 球 物理 参数 之 间 存 在 内 
在 的 制约 关系 ， 可 以 利用 卫星 参数 之 间 的 相互 制约 
关系 减少 一 个 未 知 数 we20。 从 图 2 可 以 看 出 ,不同 


B4(T;) = (1 — c; (0))(1 (à = 15(8)) B4 (Ta) + Ba(Ts)Ta(O0)ea (3) 
(BN “未知 数 ) kno CRIN) RM) 


BS RI EH IY ZY SEK), AO. NE. BA 
TE LB oy I , ERKI KIZEE REAN E AS 
化 ， 从 而 改变 发 射 率 ， 发 射 率 变化 影响 地 表 的 辐射 
效率 ， 而 地 表 温 度 变 化 决定 土壤 水 分 的 蒸发 速度 ， 
从 而 影响 与 近 地 表 空气 的 能 量 交换 和 改变 近 地 表 空 
气温 度 以 及 大 气 剖 面 温度 。 理 论 上 讲 ， 土 壤 水 分 的 
变化 也 会 影响 热 红 外 波段 发 射 率 的 变化 ， 但 传统 算 
法 中 没有 考虑 到 这 一 点 ， 一 般 认为 热 红 外 波段 的 发 
射 率 是 相对 不 变 的 ， 其 实 也 随 土壤 和 植被 含水 量 的 
变化 而 变化 。 因 此 在 降雨 比较 多 的 地 方 ， 热 红外 遥 
感 传统 方法 反 演 地 表 温 度 的 精度 波动 相对 较 大 ， 但 
AI 方法 可 以 克服 这 个 困难 。 近 地 表 空 气温 度 的 变化 
影响 大 气 剖 面 ， 从 而 影响 大 气 平均 作用 温度 。 在 地 
表 热 辐射 经 过 大 气 时 ， 被 大 气 水 汽 吸 收 ， 然 后 到 达 
Tea AK, Me EBB] AO 
单个 方程 中 有 4 个 未 知 数 。 如 果 没 有 先 验 知识 ， 必 
须 至 少 有 4 个 热 红 外 窗口 波段 从 理论 上 构建 反 演 方 
fü; 如 果 有 大 和气 水 汽 作为 先 验 知识 ， 也 可 以 只 需要 
3 个 热 红 外 波段 。 以 往 大 部 分 研究 主要 是 针对 单个 
参数 进行 反 演 ， 由 于 不 同 参数 之 间 相 互 纠 缠 ， 利 用 
AI 方法 进行 联合 反 演 和 交叉 迭代 可 以 提高 反 演 精 
度 “"。 在 进行 物理 逮 辑 推理 后 确定 物理 方法 的 理 
论 构建 后 ， 在 实际 应 用 研究 中 发 现 ， 只 利用 物理 模 
型 建立 的 方法 精度 不 是 十 分 可 靠 ， 因 为 物理 模型 对 
真实 世界 进行 了 简化 ， 而 且 很 多 地 学 模型 也 只 适合 
部 分 情况 "3。 因 此 ,为 了 尽 可 能 地 满足 更 多 的 情 
况 ， 还 要 基于 物理 逻辑 推理 进一步 建立 泛 化 的 统计 
方法 ， 这样 才 能 保证 方法 的 完整 性 和 精度 的 可 
靠 性 T, 

MODIS 卫星 传感器 是 目前 热 红 外 波段 数量 和 
质量 最 好 的 传 感 费 之 一 ， 本 研究 将 针对 MODIS 的 
热 红 外 波段 进行 模拟 和 论证 分 析 。 选 定 MODIS 波 
段 27/28/29/31/32/33 (6.5~13.5 um) ， 其 中 MODIS 
波段 27、28 和 29 在 是 热 红 外 波段 的 水 汽 吸 收 区 ， 
波段 31 和 32 在 热 红 外 波段 的 窗口 区 ，33 波段 在 热 
红外 波段 窗口 区 的 边缘 上 。 将 利用 JPL (http://spe- 
clib.jplnasa.gov) 测量 的 地 物 波 普 曲线 在 MODIS 波 
Et 27/28/29/31/32/33 的 发 射 率 作为 MODTRAN4 的 
输入 和 参数， 地 表 温 度 变化 范围 为 273~325 K， 近 地 
表 空 气温 度 变化 范围 为 273~320 K， 大 气 水 汽 含量 
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变化 范围 0.1~4.0 gem, BEIM fg Jy O~45°, 4 
模拟 观测 角度 过 大 时 ， 由 于 和 斜 路 径 上 的 水 汽 过 大 
导致 透 过 率 过 低 ， 因 此 对 这 部 分 模拟 数据 进行 了 
去 除 。 将 模拟 数据 随机 分 成 两 部 分 : 训练 数据 
45,650 组， 测试 数据 16,550 组 。 根 据 传感器 上 获得 
不 同 参数 信息 量 的 多 少 ， 将 不 同 的 参数 反 演 划分 成 
三 组 ， 分 别 对 地 球 物 理 参数 理论 和 判别 条 件 进 行 论 
证 分 析 。 第 一 组 反 演 地 表 温 度 和 发 射 率 ， 主 要 分 析 
和 论证 输入 变量 和 输出 变量 之 间 的 方程 组 个 数 和 因 
果 关 系 对 反 演 精度 影响 ; 第 二 组 反 演 近 地 表 空 气温 
度 ， 主 要 分 析 和 论证 利用 输入 和 输出 变量 之 间 强 相 
关 性 变量 作为 先 验 知识 ， 提 高 弱 关 系 反 演变 量 的 精 
度 ; 第 三 组 反 演 大 气 水 汽 含量 ， 主 要 分 析 论 证 当 输 
入 和 输出 变量 之 间 存 在 强 相 关 性 时 ， 添 加 先 验 知识 
可 以 提高 反 演 稳定 性 ， 但 对 精度 提高 不 明显 。 


收 波段 ， 增 加 热 红 外 水 汽 吸 收 波段 可 以 提高 反 演 精 
度 。 当 只 用 两 个 热 红 外 波段 (波段 31 和 32) 反 演 
时 ， 反 演 误差 接近 2K， 因 此 传统 劈 窗 算法 必须 获 
得 大 气 水 汽 和 地 表 发 射 率 作为 先 验 知识 。 而 反 演 精 
度 要 达到 1 开 以 内 ， 至 少 需要 4 个 热 红 外 波段 ， 或 
者 2 个 热 红外 波段 加 高 精度 的 地 表 发 射 率 和 水 汽 作 
为 先 验 知识 ， 也 就 证 明了 基于 AI 的 热 红 外 遥感 地 
表 温 度 高 精度 反 演 需要 理论 上 满足 方程 闭合 的 
条 件 。 

表 2 波段 28-29-31-32 组 合 时 地 表 温 度 的 反 演 误差 

Table 2 Retrieval errors of surface temperature for band 


28-29-31-32 combination 


隐 含 节点 


600 700 800 900 


NI Y & 


M/K SD/K R M/KSD/K R M/KSD/K R M/KSD/K R 


3.2 地 表 温 度 和 发 射 率 反 演 
在 3.1 节 的 物理 逻辑 推理 中 可 知 ， 反 演 地 表 温 


7 
8 
9 
10 


0.76 0.68 0.987 0.89 1.13 0.984 0.77 0.71 0.988 0.66 0.68 0.987 
0.97 2.25 0.956 0.82 0.84 0.986 0.53 0.58 0.991 0.71 0.72 0.985 
0.79 0.78 0.986 0.82 0.79 0.986 0.45 0.53 0.998 0.78 0.73 0.981 
0.83 0.89 0.986 0.71 0.78 0.987 1.03 2.63 0.928 1.14 1.56 0.963 


度 至 少 需要 4 个 热 红 外 波段 组 成 辐射 传输 方程 组 ， 
意味 着 当 输 入 热 红 外 波段 少 于 4 个 时 ， 反 演 精 度 不 
高 。 表 1 是 MODIS 波 段 29、31 和 32 组 合 在 给 定 条 
件 下 的 理论 反 演 精度 。 从 表 1 中 可 以 看 出 ， 当 隐 含 
层 是 9 层 且 每 层 隐 含 节 点 都 为 700 时 ， 平 均 理 论 精 
度 最 高 为 1.13 K， 标 准 差 为 1.17 K， 相 关系 数 为 
0.988。 
表 1 波段 29-31-32 组 合 时 地 表 温度 反 演 误差 


Table 1 Retrieval errors of surface temperature for band 


29-31-32 combination 


隐 含 节点 


np BB 


600 700 800 900 


là MKSDK R M/KSD/K R M/KSD/K R M/KSD/K R 


7 1.51 1.41 0.912 1.32 1.25 0.985 1.26 1.84 0.95 1.36 1.23 0.961 
8 1.42 1.32 0.925 1.22 1.21 0.985 1.24 1.81 0.951 1.61 1.56 0.960 
9 1.33 1.25 0.932 1.13 1.17 0.988 1.21 1.32 0.963 1.37 1.36 0.962 
10 1.27 2.04 0.912 1.23 2.48 0.910 2.46 3.11 0.879 2.34 3.77 0.896 


TE: M 为 平均 误差 ;SD 为 标准 偏差 ; R 为 相关 系数 
表 2 是 MODIS 波段 28-29-31-32 组合 反 演 的 误 
差 信息 表 。 从 表 2 中 可 以 看 出 ， 当 隐 含 层 是 9 并 且 


表 3 是 MODIS 波段 27-28-29-31-32 波段 组 合 时 
反 演 的 地 表 温 度 误差 信息 表 ， 从 表 3 中 可 以 看 出 当 
隐 含 层 数 9、 每 层 隐 含 节点 是 900 时 ， 反 演 精 度 最 
高 。 平 均 最 高 精度 是 0.44 K, 标准 差 是 0.52 K， 相 
关系 数 是 0.999。MODIS 27 波 段 也 是 水 汽 波段 ， 增 
加 水 汽 吸收 波段 能 够 提高 水 汽 的 精度 ， 理 论 上 可 以 
提高 整体 反 演 精度 。 这 里 精度 提高 不 明显 ， 主 要 原 
因 是 波段 27 的 地 面 辐射 几乎 很 难 穿 透 达 到 传感器 
也 就 说 透 过 率 非常 低 ， 星 上 亮 温 主要 包含 的 是 高 空 
大 气 水 汽 信息 ， 因 此 对 地 表 温 度 反 演 精度 贡献 不 是 
非常 大 ， 但 增加 了 反 滨 的 稳定 性 。 

表 4 是 MODIS 波段 27-28-29-31-32-33 作为 输 
人 参数 时 的 反 演 误差 信息 表 ， 当 隐 含 层 是 8 并 且 每 
表 3 波段 27-28-29-31-32 组 合 时 地 表 温 度 的 反 演 误差 


Table 3 Retrieval errors of surface temperature for band 27-28- 


29-31-32 combination 


隐 含 节点 


600 700 800 900 


M/K SD/K R M/KSD/K R M/KSDK R M/KSD/K R 


每 层 节 点 为 800 时 反 演 精度 最 高 。 最 高 平均 精度 
为 0.45 K, PEZH 0.53 K 和 相关 系数 为 0.998。 
从 表 1 和 表 2 可 以 看 出 ， 当 用 4 个 热 红 外 波段 时 ， 
平均 精度 提高 了 0.68 Ko MODIS 波段 28 是 水 汽 吸 


o6 o x WI] 2» 88 


0.64 0.68 0.993 0.58 0.59 0.995 0.61 0.69 0.995 0.76 0.88 0.979 
0.61 0.67 0.995 0.62 0.93 0.993 0.66 0.85 0.978 0.48 0.54 0.998 
0.62 0.68 0.994 0.65 1.02 0.991 0.73 0.97 0.965 0.44 0.52 0.999 
0.65 0.88 0.978 0.64 1.12 0.99 0.48 0.61 0.998 0.56 0.89 0.997 
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层 隐 含 节点 为 900 时 精度 最 高 。 平 均 精 度 最 高 为 
0.51 K， 标 准 差 为 0.55 K， 相 关系 数 是 0.997。 从 
表 3 和 表 4 可 以 看 出 ， 当 有 汽 波段 时 ， 添 加 波段 33 
时 并 没有 提高 反 演 精度 。 主 要 原因 是 热 红 外 波段 33 
受 CO, 的 影响 比较 大 ， 不 是 非常 适合 于 反 演 地 表 温 
度 ， 因 此 增加 热 红 外 波段 不 一 定 会 增加 反 演 精度 ， 
有 时 反而 会 降低 反 演 精度 。 表 3 分 析 说 明 当 增加 的 
波段 与 输出 参数 之 间 存 在 强 相 关 性 时 ， 可 以 提高 精 
度 ; 表 4 说 明 当 增加 的 波段 与 输出 参数 之 间 是 弱 关 
系 或 者 是 增加 了 噪声 时 ， 会 降低 反 演 精度 。 

同时 也 可 以 将 发 射 率 反 演 计算 出 来 ， 这 里 主要 
列 出 了 波段 27-28-29-31-32 组 合 时 波段 31 和 32 的 发 
射 率 误 差 ， 如 表 5 和 表 6 所 示 ， 反 演 误 差 都 在 0.01 


表 4 波段 27-28-29-31-32-33 组 合 时 地 表 温 度 的 反 演 
误差 
Table 4 Retrieval errors of surface temperature for band 27-28- 
29-3 1-32-33 combination 


隐 含 节点 
600 700 800 900 
Je M/K SD/K R M/KSD/K R M/KSD/K R M/KSD/K R 


zp 8E 


7 0.68 0.69 0.991 0.65 0.67 0.991 0.65 0.68 0.996 0.63 1.21 0.957 
8 0.62 0.65 0.992 0.63 0.88 0.992 0.65 0.70 0.995 0.51 0.55 0.997 
9 0.64 0.66 0.991 0.68 0.68 0.995 0.69 0.79 0.992 0.54 0.56 0.996 
10 0.93 1.65 0.912 0.69 0.72 0.994 0.51 0.55 0.998 0.77 1.51 0.935 


以 下 。 发 射 率 在 热 红 外 波段 相对 被 动 微 波 受 土壤 水 
分 变化 的 影响 较 小 ， 能 够 用 来 表征 地 表 类 型 。 在 微 
波 波段 ， 发 射 率 受 土壤 水 分 的 变化 影响 很 大 ， 地 表 
温度 和 土壤 水 分 通过 发 射 率 纠缠 在 一 起 " 。 


表 5 波段 27-28-29-31-32 组 合 时 波段 31 的 发 射 率 误差 
Table 5 Retrieval emissivity errors in band 31 for band 27-28-29-31-32 combination 


CERIS 
隐 含 层 600 700 800 900 
M/K SD/K R M/K SD/K R M/K SD/K R M/K SD/K R 
0.006 0.007 0.971 0.006 0.007 0.976 0.008 0.086 0.926 0.005 0.006 0.981 
8 0.007 0.008 0.965 0.005 0.006 0.980 0.005 0.007 0.979 0.007 0.023 0.953 
9 0.005 0.007 0.972 0.005 0.007 0.978 0.008 0.021 0.944 0.005 0.006 0.983 
10 0.005 0.006 0.976 0.004 0.007 0.981 0.004 0.005 0.991 0.007 0.007 0.961 
HO 波段 27-28-29-31-32 组 合 时 波段 32 的 发 射 率 误差 
Table 6 Retrieval emissivity errors in band 32 for band 27-28-29-31-32 combination 
隐 含 节点 
隐 含 层 600 700 800 900 
M/K SD/K R M/K SD/K R M/K SD/K R M/K SD/K R 
0.005 0.006 0.972 0.005 0.006 0.977 0.005 0.007 0.976 0.006 0.008 0.962 
8 0.005 0.006 0.975 0.004 0.005 0.986 0.004 0.005 0.986 0.004 0.004 0.992 
0.004 0.005 0.979 0.005 0.013 0.961 0.008 0.014 0.948 0.005 0.005 0.981 
10 0.004 0.006 0.978 0.004 0.006 0.977 0.007 0.006 0.956 0.005 0.006 0.977 


3.9 近 地 表 空气 温度 反 演 


近 地 表 空气 温度 是 天 气 预 报 和 农业 干旱 等 模型 
的 关键 参数 ， 以 往 主 要 是 通过 气象 站 点 插值 获取 。 
也 有 一 些 研究 通过 统计 方法 或 者 神经 网 络 从 遥感 数 
据 中 反 演 近 地 表 空气 温度 ， 但 通常 缺乏 可 解释 和 物 
理 意义 ， 而 且 精 度 有 限 。 本 研究 采用 深度 学 习 和 基 
于 物理 逻辑 推理 反 演 近 地 表 空气 温度 ， 其 推理 和 分 
析 过 程 参 考 文献 [9] 和 [15]。 在 热 红 外 遥感 窗口 
波段 ， 热 红外 传感器 获取 的 主要 信息 来 自 地 表 ， 近 
地 表 空 气温 度 的 信息 量 相 对 较 少 ， 因 此 反 演 精度 受 


限 。 表 7 是 MODIS 波段 27-28-29-31-32 组 合 时 直接 
反 演 近 地 表 空气 温度 的 误差 信息 表 ， 当 隐 含 层 是 10 
和 每 层 隐 含 节 点 是 700 时 精度 最 高 。 理 论 平 均 精 度 
最 高 为 1.42 K， 标 准 差 为 1.46 K， 相 关系 数 0.975。 

K 8 是 波段 27-28-29-31-32+LST+LSE31+ 
LSE32 组 合 时 近 地 表 空气 温度 的 反 演 误 差 (其 中 ， 
LSE31 为 波段 31 发 射 率 ，LSE32 为 波段 32 发 射 率 )， 
最 高 精度 是 在 隐 含 层 为 10 和 每 层 隐 舍 节点 为 800 
时 ， 最 高 平均 精度 为 0.81 K， 标 准 偏差 为 0.91 K, 
相关 系数 为 0.984。 当 利用 地 表 温 度 和 发 射 率 作 为 
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#7 波段 27-28-29-31-32 组 合 时 近 地 表 空气 温度 反 演 
误差 
Table 7 Retrieval errors of near-surface air temperature for 


band 27-28-29-3 1-32 combination 


表 9 波段 27-28-29-31-32 组 合 时 大 气 水 汽 含 量 反 演 误 差 


Table 9 Retrieval errors of atmospheric water vapor content for 


band 27-28-29-31-32 combination 


隐 隐 含 节点 
含 600 700 800 900 
Ez M/KSDK R M/KSDK R M/KSD/K R MKSDK R 


600 700 800 900 


Ml mb 8 


M/K SD/K R M/K SDK R M/KSD/K R M/KSD/K R 


7 1.64 1.55 0.958 1.52 1.55 0.962 1.58 1.66 0.963 1.68 1.78 0.966 
8 1.58 1.64 0.953 1.56 1.58 0.958 1.56 1.65 0.964 1.60 1.72 0.958 
9 1.44 1.51 0.959 1.49 1.54 0.965 1.56 1.63 0.968 1.69 1.77 0.959 
10 1.47 1.49 0.961 1.42 1.46 0.975 1.46 1.52 0.967 1.64 1.69 0.969 


先 验 知识 时 ， 近 地 表 空 气温 度 的 精度 明显 得 到 提 
高 ， 而 且 反 演 精度 更 加 稳定 。 这 主要 是 因为 利用 地 
表 温 度 作 为 先 验 知识 放大 了 近 地 表 空气 温度 的 信 
号 ， 而 且 发 射 率 提 供 了 地 表 类 型 ， 因 此 反 演 精度 大 
幅 提 高 ， 并 且 算 法 更 具 移 植 性 ”“。 因 此 对 于 输入 变 
量 与 输出 变量 之 间 是 弱 相 关 的 参数 反 演 ， 加 入 强 相 
关 性 的 变量 作为 先 验 知识 可 以 提高 弱 相 关 参 数 的 反 
演 精度 。 

表 8 波段 27-28-29-31-32+LST+LSE31+LSE32 组 合 时 

近 地 表 空气 温度 反 演 误差 


Table 8 Retrieval errors of near-surface air temperature for 


band 27-28-29-3 1-32+LST+LSE31+LSE32 combination 


7 0.18 0.19 0.960 0.15 0.19 0.971 0.18 0.21 0.977 0.15 0.17 0.979 
8 0.12 0.15 0.975 0.23 0.27 0.972 0.11 0.13 0.979 0.13 0.15 0.980 
9 0.17 0.22 0.963 0.18 0.23 0.973 0.09 0.11 0.989 0.14 0.16 0.976 


10 0.87 0.93 0.875 0.35 0.41 0.938 0.13 0.15 0.983 0.40 0.55 0.923 


隐 隐 含 节点 
Tr 600 700 800 900 
Ez M/KSD/K R M/KSD/K R M/KSDK R M/KSD/K R 


7 1.14 1.36 0.976 1.19 1.37 0.975 1.14 1.38 0.981 1.16 1.36 0.978 
8 1.25 1.34 0.968 1.18 1.38 0.975 1.19 1.41 0.978 1.03 1.16 0.979 
9 1.22 1.32 0.971 1.26 1.41 0.964 0.93 1.05 0.980 1.19 1.37 0.975 
10 1.19 1.29 0.976 1.14 1.36 0.977 0.81 0.91 0.984 1.22 1.34 0.973 


3.4 大 气 水 汽 含 量 反 演 

大 气 水 汽 含量 是 天 气 预报 和 农业 干旱 监测 等 
模型 的 重要 参数 ， 在 上 面 分 析 的 基础 上 ， 表 9 是 波 
段 27-28-29-31-32 组 合 时 的 大 气 水 汽 含量 反 演 误 
差 。 最 高 精度 是 在 隐 含 层 数 为 9 和 每 层 隐 含 节点 为 
800 时 ， 最 高 平均 精度 为 0.09 g/cm*， 标准 偏差 为 
0.11 g/cm ， 相 关系 数 为 0.989。 

表 10 £2 YK Bt 27-28-29-31-32+ LST+LSE 组 合 时 
的 大 气 水 汽 含量 反 演 误差 ， 最 高 精度 是 在 隐 含 层 
为 9 和 每 层 隐 含 节点 为 900 时 ， 最 高 平均 精度 为 
0.08 g/cm*， 标 准 偏 差 为 0.09 gcm， 相 关系 数 为 


0.992。 对 比 表 9 和 表 10， 当 使 用 LST 和 LSE 作 为 先 
验 知识 时 ， 精 度 略 有 提高 和 误差 相对 比较 稳定 。 如 
果 当 先 验 知识 LST 和 LSE 具 有 一 定 误 差 时 ， 可 能 精 
度 不 会 有 太 大 的 改善 。 因 此 ， 当 输入 参数 中 有 2 个 
输入 波段 对 输出 参数 水 汽 比较 敏感 时 ， 可 以 不 加 入 


先 验 知识 。 

表 10 波段 27-28-29-31-32+LST+LSE 组 合 时 大 气 水 汽 
含量 的 反 演 误差 

Table 10 Retrieval errors of atmospheric water vapor content 


for band 27-28-29-31-32+LST+LSE combination 


隐 含 节点 
600 700 800 900 
Je M/K SDK R MKSD/K R M/KSD/K R M/KSD/K R 
7 0.15 0.17 0.976 0.18 0.21 0.971 0.17 0.18 0.975 0.14 0.16 0.976 
8 0.15 0.16 0.976 0.21 0.23 0.969 0.09 0.11 0.991 0.15 0.17 0.977 
9 0.14 0.16 0.977 0.23 0.25 0.967 0.11 0.15 0.979 0.08 0.09 0.992 


10 0.17 0.19 0.975 0.15 0.19 0.971 0.09 0.11 0.990 0.11 0.12 0.989 


n BS 


4 讨 论 


人 类 认识 世界 通过 观察 (取样 )， 然 后 基于 观 
察 和 思考 用 大 脑 抽 象 构造 模式 ， 最 大 程度 地 形成 知 
识 以 达到 快速 认识 世界 或 者 目标 的 目的 ， 本 质 上 可 
以 归结 为 定义 规则 认识 世界 。 目 标 只 有 一 个 一 一 最 
大 程度 地 提高 识别 〈 认 知 ) 精度 ， 因 此 不 必 过 于 考 
虑 具体 方法 的 形式 〈 统 计 方 法 或 物理 方法 或 AI 方 
法 )， 这 些 都 是 认 知 的 手段 。 根 据 作 者 团队 多 年 的 
研究 经 验 ， 统 计 方法 、 物 理 方法 和 机 器 学 习 方 法 本 
质 上 是 一 致 的 ， 是 人 们 根据 不 同 的 情况 定义 的 认识 
事物 规律 的 一 种 方法 或 手段 ， 特 别 是 目前 用 得 比较 
多 的 深度 学 习 神 经 网 络 算法 其 实 是 一 种 更 高 级 的 统 
计 优 化 计算 方法 ， 这 种 方法 的 精度 和 普 适 范围 取决 
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于 构建 样本 库 (训练 和 测试 数据 的 代表 性 和 精度 ) 
时 遵循 的 规则 是 否 更 符合 真实 情况 。 

传统 遥感 地 球 物理 参数 反 演 方法 主要 适合 于 地 
势 平 坦 ， 地 表 类 型 单一 的 纯净 像 元 。 如 图 3 所 示 ， 
对 于 地 球 物理 遥感 参数 反 演 , 传统 的 物理 方法 无 
法 解决 像 元 相对 入 射 角 问 题 ， 也 很 难 解 决 混合 像 
元 问题 ， 从 而 导致 传统 方法 反 演 精度 具有 瓶颈 ， 公 
X (4) 中 参数 的 物理 意义 具体 请 参考 文献 [9] 和 
[15]。AI 方 法 可 以 将 混合 像 元 当成 一 个 新 的 等 效 地 
物 类 型 ， 在 解决 复杂 像 元 反 演 问题 上 具有 比 传统 物 
理 方法 非常 强 的 优势 。 从 能 量 信息 流传 播 的 角度 
看 ， 深 度 学 习 神 经 网 络 通过 神经 元 控制 传递 能 量 的 
方式 比 统计 和 物理 方法 更 接近 真实 情况 ， 甚 至 更 具 
有 物理 意义 ， 这 需要 人 们 提高 或 改变 认 知 思维 ， 并 
从 信息 传播 的 角度 解读 天。 因此 只 要 深度 学 习 耦 
合 物理 和 统计 方法 的 模式 满足 范式 条 件 ， 理 论 上 可 
以 无 限 逼 近 真 实情 况 。 


100 m 
大 下 降 
气 对 流 周 0.6 K 


(3 zi 
By(T;) = t,(8)[|5(8)BCT,) + (1 — &(0))f(Taj)| + f(Tad) (4) 
图 3 地 球 物 理 参 数 逐 感 反 演 混合 像 元 问题 

Fig. 3 Mixed pixel problem in remote sensing retrieval of geo- 


physical parameters 


另外 ， 从 某 个 角度 上 讲 ，AI 是 一 种 高 级 的 优化 
计算 方法 ， 与 物理 方法 和 统计 方法 在 本 质 上 是 一 样 
的 ， 都 是 识别 目标 的 手段 。 只 不 过 机 器 学 习 方 法 在 
硬件 技术 高 速 发 展 的 辅助 下 ， 显 示 出 更 加 强大 的 识 
别 目标 的 迭代 优化 计算 能 力 ， 深 度 学 习 理 论 上 可 以 
将 很 多 组 不 同 物理 方法 和 统计 方法 耦合 在 一 起 。 在 
满足 反 演 范式 条 件 下 ， 深 度 学 习 的 反 演 和 适应 范围 
理论 上 远 远 胜 过 一 个 或 者 几 个 物理 和 统计 传统 方法 
的 组 合 。 作 者 团队 在 2006 年 和 2007 年 将 深度 学 习 
(动态 学 习 神 经 网 络 ) 应 用 于 地 表 温 度 和 发 射 率 反 
演 "” ,但 当时 由 于 计算 机 内 存 不 够 ， 隐 含 层 一 般 


只 设置 了 3 层 ， 而 现在 算法 的 隐 含 层 和 内 存 翻 到 了 
几 信 甚至 几 十 倍 上 百倍 "“ “， 对 计算 机 要 求 更 高 。 
由 于 计算 机 技术 的 进步 ， 目 前 地 球 物 理 参数 遥感 反 
演 的 精度 也 大 幅 提 升 。 在 满足 范式 条 件 下 ， 基 于 AI 
的 深度 学 习 耦 合 物理 和 统计 方法 范式 是 目前 利用 遥 
感 数据 反 演 地 球 物理 参数 特别 是 农业 气象 遥感 关键 
参数 的 最 佳 方法 。 

本 文 主要 讨论 的 是 地 球 物理 参数 反 演 范式 理论 
和 判定 条 件 ， 其 对 于 遥感 分 类 或 者 目标 识别 也 同样 
i&H ^", BRERA (Convolutional Neu- 
ral Network, CNN) 被 认为 是 较 好 的 分 类 和 目标 识 
别 方法 之 一 ， 其 实 CNN 是 在 传统 神经 网 络 的 基础 
上 进一步 利用 不 同 的 卷 积 核 从 不 同 的 尺度 上 提取 目 
标的 信息 ， 目 的 就 是 通过 不 同 维度 的 信息 来 唯一 确 
定 目标 ， 从 而 达到 提高 分 类 (目标 识别 ) 的 精 
度 将。 理论 上 讲 ， 如 果 这 些 从 不 同 维度 (不 同 卷 
积 核 ) 上 提取 的 信息 也 可 以 构建 数学 方程 唯一 确定 
输出 参数 ， 也 可 以 形成 范式 ,否则 同样 也 存在 不 确 
定性 。 


5 结 论 


最 近 几 年 ，AI 在 许多 领域 的 突破 催生 了 很 多 新 
技术 、 新 产品 和 新 产业 ， 将 深刻 影响 和 改变 人 类 生 
产生 活 方式 和 思维 模式 ， 实 现 生产 力 的 整体 跃升 。 
为 了 促进 AI 在 地 学 和 农学 的 应 用 和 培养 交叉 学 科 
的 人 才 ， 本 研究 提出 了 基于 AI 地 球 物 理 参 数 反 演 
范式 理论 及 判定 条 件 ， 即 通过 物理 逻辑 推理 ， 目 标 
信息 (问题 ) 在 理论 或 物理 逻辑 上 可 以 用 数学 表达 
式 进行 描述 ， 同 时 可 以 用 物理 模型 模拟 获得 物理 方 
法 的 解 。 为 了 克服 物理 方法 的 缺陷 ， 基 于 物理 逻辑 
推理 构建 泛 化 的 统计 方法 ， 并 且 通 过 多 源 数据 获得 
统计 方法 的 解 。 在 此 基础 上 ， 物 理 方 法 和 统计 方法 
的 解构 成 深度 学 习 的 训练 和 测试 数据 ， 从 而 达到 深 
度 学 习 炮 合 物理 和 统计 方法 的 目的 。 深 度 学 习 通 过 
大 数据 学 习 和 优化 计算 将 物理 方法 和 统计 方法 相 结 
合 ， 从 而 保证 了 提出 的 范式 不 仅 有 物理 意义 和 可 解 
释 性 ， 而 且 还 具备 通用 性 。 

AI 方法 反 演 地 球 物理 参数 是 否 形 成 范式 ， 即 形 
成 具有 物理 机 理 和 可 解释 的 次 度 学 习 范 式 基 本 条 件 
有 两 个 : (1) 输入 与 输出 变量 (参数) 之 间 必 须 存 
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在 因果 关系 ; (2) 输入 和 输出 变量 (参数 ) 之 间 理 
论 上 可 以 构建 闭合 的 方程 组 (未知 数 个 数 少 于 或 等 
于 方程 组 个 数 )， 也 就 是 说 输出 参数 可 以 被 输入 变 
量 唯 一 确定 。 如 果 使 用 深度 学 习 没有 证 明 或 不 能 证 
明 存 在 上 述 两 个 条 件 ， 一 般 被 认定 为 “黑箱 ”。 如 
果 能 证 明 : 输入 和 输出 参数 之 间 肯 定 不 存在 闭合 的 
关系 (未知数 大 于 方程 组 个 数 )， 那 么 这 样 的 应 用 
可 能 主要 适合 于 局 部 区 域 ， 不 能 通用 ， 很 难 具备 移 
植 性 ， 不 能 称 为 范式 。 

AT 方法 的 应 用 不 能 脱离 具体 的 问题 ， 在 解决 某 
个 应 用 问题 时 ， 要 先 对 问题 对 象 进行 物理 逻辑 推 
理 ， 从 理论 上 构造 方程 ， 确 定 是 输入 变量 和 输出 变 
量 之 间 的 因果 关系 ; 然后 从 方程 里 面 的 未 知 数 个 数 
确定 方程 组 的 个 数 ， 也 就 是 确定 几 个 输入 参数 ， 最 
后 使 用 深度 学 习 进 行 优 化 求解 。 进 行 物 理 逻 辑 推理 
是 使 得 深度 学 习 具 有 物理 意义 和 可 解释 的 关键 。 利 
用 MODIS 遥感 数据 反 演 地 表 温 度 、 发 射 率 、 近 地 
表 空 气温 度 和 大 气 水 汽 含量 对 范式 理论 进行 了 论 
WE: 当 输 出 参数 (LST 和 LSE) 与 输入 变量 (BT) 
之 间 具 有 强 相 关 性 时 ， 利 用 深度 学 习 耦 合 物理 和 统 
计 方 法 可 以 获得 非常 高 的 精度 ; 当 输 出 参数 
(NSAT) 与 输入 变量 (BT) 之 间 是 弱 相 关 性 时 ， 
加 入 先 验 知识 (LST 和 LSE) 可 以 提高 输出 参数 
(NSAT) 反 演 精度 和 稳定 性 ; 部 分 强 相 关 性 时 
(WVC 与 BT,)， 加 入 先 验 知识 (LST 和 LSE) 虽然 
能 稍微 提高 精度 和 稳定 性 ， 但 先 验 知识 (LST 和 
LSE) 的 误差 可 能 带 来 不 确定 性 ， 因 此 也 可 以 不 加 
入 先 验 知 识 。 通 过 对 MODIS 传感器 热 红 外 波段 地 
球 物 理 参 数 反 演 分 析 可 知 ， 波 段 27、28、29、31 
更 适合 于 反 演 大 气 水 汽 含量 ， 波 段 28、29、31 和 
32 更 适合 于 反 演 地 表 温 度 和 发 射 率 及 近 地 表 空气 温 
度 。 如 果 要 获得 4 个 参数 的 最 高 精度 ， 则 建议 仪器 
设计 5 个 波段 (27、28、29、31、32) 最 合适 。 如 
果 只 设计 4 个 热 红 波段 ， 可 以 优先 考虑 27、28、 
31、32。 本 研究 分 析 结 果 表 明 地 球 物 理 参 数 反 演 范 
式 理论 和 判定 条 件 可 靠 ， 该 理论 和 判定 条 件 的 提出 
对 利用 AI 方法 反 演 地 表 物 理 参 数 具 有 重大 意义 。 

地 球 物理 参数 反 演 范式 理论 和 判定 条 件 也 适用 
于 遥感 分 类 等 目标 识别 ， 只 是 需要 换 一 个 角度 解 
读 ， 比 如 不 同 卷 积 核 提 取 的 特征 信息 要 能 唯一 确定 


目标 。 在 满足 范式 条 件 下 ， 基 于 AI 反 演 地 球 物理 
参数 是 最 佳 选 择 ， 因 此 建议 科技 部 和 国家 自然 基金 
委 给 予 支 持 ， 打 造 中 国 基 于 AI 的 地 球 物理 参数 反 
演 范 式 的 “ChatGPT”。 


利益 冲突 声明 : 本 研究 不 存在 研究 者 以 及 与 公开 
研究 成 果 有 关 的 利益 冲突 。 
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[Objective] Deep learning is one of the most important technologies in the field of artificial intelligence, which has sparked a re- 
search boom in academic and engineering applications. It also shows strong application potential in remote sensing retrieval of geo- 


physical parameters. The cross-disciplinary research is just beginning, and most deep learning applications in geosciences are still 
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"black boxes", with most applications lacking physical significance, interpretability, and universality. In order to promote the applica- 
tion of artificial intelligence in geosciences and agriculture and cultivate interdisciplinary talents, a paradigm theory for geophysical 
parameter retrieval based on artificial intelligence coupled physics and statistical methods was proposed in this research. 

[Methods] The construction of the retrieval paradigm theory for geophysical parameters mainly included three parts: Firstly, physical 
logic deduction was performed based on the physical energy balance equation, and the inversion equation system was constructed the- 
oretically which eliminated the ill conditioned problem of insufficient equations. Then, a fuzzy statistical method was constructed 
based on physical deduction. Representative solutions of physical methods were obtained through physical model simulation, and oth- 
er representative solutions as the training and testing database for deep learning were obtained using multi-source data. Finally, deep 
learning achieved the goal of coupling physical and statistical methods through the use of representative solutions from physical and 
statistical methods as training and testing databases. Deep learning training and testing were aimed at obtaining curves of solutions 
from physical and statistical methods, thereby making deep learning physically meaningful and interpretable. 

[Results and Discussions] The conditions for determining the formation of a universal and physically interpretable paradigm were: 
(1) There must be a causal relationship between input and output variables (parameters); (2) In theory, a closed system of equations 
(with unknowns less than or equal to the number of equations) can be constructed between input and output variables (parameters), 
which means that the output parameters can be uniquely determined by the input parameters. If there is a strong causal relationship be- 
tween input parameters (variables) and output parameters (variables), deep learning can be directly used for inversion. If there is a 
weak correlation between the input and output parameters, prior knowledge needs to be added to improve the inversion accuracy of 
the output parameters. The MODIS thermal infrared remote sensing data were used to retrieve land surface temperature, emissivity, 
near surface air temperature and atmospheric water vapor content as a case to prove the theory. When there was strong correlation be- 
tween output parameters (LST and LSE) and input variables (BTi), using deep learning coupled with physical and statistical methods 
could obtain very high accuracy. When there was a weak correlation between the output parameter (NSAT) and the input variable 
(BTi), adding prior knowledge (LST and LSE) could improve the inversion accuracy and stability of the output parameter (NSAT). 
When there was partial strong correlation (WVC and BTi), adding prior knowledge (LST and LSE) could slightly improve accuracy 
and stability, but the error of prior knowledge (LST and LSE) may bring uncertainty, so prior knowledge could also be omitted. Ac- 
cording to the inversion analysis of geophysical parameters of MODIS sensor thermal infrared band, bands 27, 28, 29 and 31 were 
more suitable for inversion of atmospheric water vapor content, and bands 28, 29, 31 and 32 were more suitable for inversion of sur- 
face temperature, Emissivity and near surface air temperature. If someone want to achieve the highest accuracy of four parameters, it 
was recommended to design the instrument with five bands (27, 28, 29, 31, 32) which were most suitable. If only four thermal infrared 
bands were designed, bands 27, 28, 31, and 32 should be given priority consideration. From the results of land surface temperature, 
emissivity, near surface air temperature and atmospheric water vapor content retrieved from MODIS data using this theory, it was not 
only more accurate than traditional methods, but also could reduce some bands, reduce satellite load and improve satellite life. Espe- 
cially, this theoretical method overcomes the influence of the MODIS official algorithm (day/night algorithm) on sudden changes in 
surface types and long-term lack of continuous data, which leads to unstable accuracy of the inversion product. The analysis results 
showed that the proposed theory and conditions are feasible, and the accuracy and applicability were better than traditional methods. 
The theory and judgment conditions of geophysical parameter retrieval paradigms were also applicable for target recognition such as 
remote sensing classification, but it needed to be interpreted from a different perspective. For example, the feature information extract- 
ed by different convolutional kernels must be able to uniquely determine the target. Under satisfying with the conditions of paradigm 
theory, the inversion of geophysical parameters based on artificial intelligence is the best choice. 

[Conclusions] The geophysical parameter retrieval paradigm theory based on artificial intelligence proposed in this study can over- 
come the shortcomings of traditional retrieval methods, especially remote sensing parameter retrieval, which simplify the inversion 
process and improve the inversion accuracy. At the same time, it can optimize the design of satellite sensors. The proposal of this theo- 


ry is of milestone significance in the history of geophysical parameter retrieval. 
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